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ABSTRACT 

We show that the re-brightening in the GRB 021004 optical afterglow light 
curve can be explained within the framework of the standard fireball model. The 
optical light curve of the forward shock emission is expected to rise initially and 
to decay after the typical synchrotron frequency crosses the optical band. We 
show that such a rising phase emission was caught for GRB 021004, together with 
a reverse shock emission. With the standard values of parameters obtained in 
other afterglow observations, we can construct example cases in which theoretical 
estimates reasonably fit the broadband observations. Therefore, the early re- 
brightening might be a common feature in optical afterglows. 



1. introduction 

GRB 021004 triggered HETE II on 2002 October 4 at 12:06:13 UT. The burst lasted 
T ~ 100 seconds (Shirasaki et al. 2002) and the 7 - 400 keV fluence was ~ 3.2 x 10~ 6 ergs 
cm" 2 (Lamb et al. 2002). The prompt localization of GRB 021004 by HETE II allowed 
the follow-up of the afterglow at very early time. Fox (2002) detected an optical transient 
~ 9 mins after the trigger at the level of R ~ 15.56 mag. Torii, Kato and Yamaoka (2002) 
observed the error box of the burst ~ 3.5 mins after the trigger, yielding upper limit around 
R ~ 13.6 mag at the position of the transient discovered by Fox (2002). The spectroscopic 
observations of the optical afterglow revealed an emission line interpreted as Ly-ct at z — 
2.328 (Mirabal et al. 2002). Assuming fi = 0.3, A = 0.7 and h = 0.6, the isotropic 
gamma-ray energy is ~ 5.6 x 10 52 ergs. 

The dense sampling of the afterglow light curve at early time revealed the peculiarity, 
a major re-brightening around ~ 0.1 day after the trigger and a short time scale variability 
around ~ 1 day. The slope of the afterglow is ~ t~ ar for the earliest three observations (9-17 
min after trigger) by Fox (2002), and after the luminosity increases around ~ 0.1 days, it 
decayed with ~ t~ im as usual afterglows do. These temporal features might be modeled 
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by refreshed shocks (Zhang and Meszaros 2002). Recently, Lazzati et al. (2002) interpreted 
the features as due to enhancements in the ambient medium density. In this Letter, we 
show that the major re-brightening around ~ 0.1 day can be explained within the standard 
fireball model. We propose that the reverse shock emission dominated the optical band at 
early times, leading to the observed peculiarity. 



2. The model 

Consider a relativistic shell (fireball ejecta) with energy E, a Lorentz factor rj and a 
width A expanding into the homogeneous ISM with a particle number density n. When the 
shell sweeps a large volume of the ISM, it is decelerated and the kinetic energy is transfered to 
the ISM by shocks (e.g. Kobayashi, Piran & Sari 1999). The shocked ISM forms a relativistic 
blast wave (forward shock) and emits the internal energy via synchrotron process. 

The emission from a reverse shock was also predicted (Meszaros & Rees 1997; Sari & 
Piran 1999a). When the reverse shock crosses the shell, the forward shocked ISM and the 
reverse shocked shell carry comparable amounts of energy. However, the typical temperature 
of the shocked shell is lower since the mass density of the shell is higher. Consequently, the 
typical frequency in the shocked shell is lower. A prompt optical emission from GRB 990123 
(Akerlof et al. 1999) can be regarded as this emission (Sari & Piran 1999b; Kobayashi & 
Sari 2000; Soderberg & Ramirez-Ruiz 2002). 



2.1. Forward Shock 

Observations of optical afterglows usually start around several hours after the burst 
trigger. Since, at such a late time, the typical synchrotron frequency of the forward shock 
emission u m j is lower than optical band u R ~ 5 x 10 14 Hz, the evolution of afterglows is 
well described by a single power law, except for the jet break (Rhoads 1999; Sari, Piran 
& Halpern 1999). However, an optical light curve is expected to peak at early time when 
the typical synchrotron frequency crosses the optical band. Before the peak time t m j, the 
luminosity increases as oc t 1//2 , and reaches to the maximum flux F U:max j, and then decays as 

t 3(l-p)/4 where 

p is the index of the power law distribution of random electrons accelerated 
at shock. Using the results in Sari, Piran and Narayan (1998), we get 

v mJ {t) ~ 5.1xl0 15 (l + z)^ e V 2 e ^VV/ 2 Hz, (1) 

u cJ (t) ~ 2.7xl0 12 (l + ^)-V 2 eB 3 / 2 ^V^ 1/2 Hz, (2) 
t mJ ~ 2.9(1 + ^)V34/3 eW /3 £ i/3 v -2/3 day8j (3) 
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F„, maxJ ~ l.lxl0 2 a u (l + z)e 1 fE 52 n 1 /2 D^ 2 mJy. (4) 

where v c j is the cooling frequency, e# and e e are the fractions of the shock energy given to 
magnetic field and electrons at the shock, g — {p — 2)/(p — 1), E 52 = E/10 52 ergs, no = n/1 
proton cm -3 , z/r 5 i5 = z/r/10 15 Hz, and here td is the observer's time in unit of day, D 2 s is the 
luminosity distance in unit of 10 28 cm, a u is a correction factor to the extinction along the 
line of sight to the burst. 

Assuming that the re-brightening with the peak luminosity of ~ 1 mJy around ~ 0.1 
day is caused by this peak, we get the following formulae for e e and no as functions of e# 
and other known parameters. 

---^^^)" 1/4 (4)" I (^) 3/4 (S)" 1/4 (^) IA <■> 

Holland et al. (2002a) claim the extinction Ay = 0.26 mag, which implies a ~ 20% correction 
in the R-band. The slope ~ —1.05 of the optical afterglow at late time implies p ~ 2.4. 
However, there is still some debate about the value of p (Sako & Harrison 2002; Holland et 
al. 2002a; Pandey et al. 2002). We will discuss two cases p = 2.2 and 2.4. 

2.2. Reverse Shock 

The evolution of reverse shocks is classified into two cases (Kobayashi 2000). If the initial 
Lorentz factor of the shell 77 is larger than a critical value r] c = (3E / '327mm p c 2 Aq) 1 / 8 where 
nip is the mass of proton, the reverse shock becomes relativistic in the frame of unshocked 
shell material during crossing the shell, and drastically decelerates the shell (thick shell case). 
If rj < r] c , the reverse shock can not decelerate the shell effectively (thin shell case). According 
to the internal shock model, the initial width of the shell A is given by the intrinsic duration 
of the GRB ~ cT/(l + z) (Kobayashi, Piran & Sari 1997). 

The Lorentz factor at the shock crossing time is given given by 7 X ~ min[r],r] c ]. The shock 
crossing time is t x ~ ("f x /ri c )~ 8 / 3 T (Sari & Piran 1999a; Kobayashi 2000). 

Since at the shock crossing time, the forward and reverse shocked regions have the same 
Lorentz factor and internal energy density e, the cooling frequency of the reverse shock z/ c r 
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is equal to that of the forward shock v c j. 

^Cr(tx) ~ Vc,f(t x ). (8) 

The typical frequency of synchrotron emission is proportional to the electron's random 
Lorentz factor squared and to the magnetic field and Lorentz boost. The Lorentz boost 
and the magnetic field oc \fe are the same for the two shocked regions, while the random 
Lorentz factor is proportional to 7x in the reverse shocked region and to 7 X in the forward 
shocked region where 7x is the Lorentz factor of the shocked shell material in the frame of 
the unshocked shell material. Using a relation 7 x 7x ~ Vi the reverse shock frequency at the 
crossing time (the peak time) is given by (Sari & Piran 1999a) 

T] 2 

Vrn,r{ty.) 4 fm,/(*x)- (9) 

7x 

The peak flux at the typical frequency is proportional to the number of electrons and to the 
magnetic field and the Lorentz boost. From the energy conservation, the mass of the shell 
is larger by a factor oi^/r] at the crossing time than that of the ISM swept by the forward 
shock. Since the number of electrons is proportional to the mass, we get 

7x 

Even though the hydrodynamic evolution of "thin" and "thick" shells are very different, 
the time dependences of the emission are similar (Kobayashi & Sari 2000). If the optical 
band vr is below the typical synchrotron frequency u mj1 ., the luminosity decays as ~ t~ ' 5 . 
If vr > v mi r, it decreases as ~ t~ 2 (Kobayashi 2000). 

First we assume that the optical band is below the typical frequency of the reverse 
shock emission at the shock crossing time v m! r(t x ) > v r ~ ^m,/(4i,/)- Using u m j oc t~ 3//2 
(Sari, Piran & Narayan 1998) and eq. (9), we get rj > {T/t mJ ) 3 / A r] 2 c ~ 0.03^ > 1000 for 
the standard values of parameters assumed in eq. (7). Then, GRB021004 should be a thick 
shell case. However, in a thick shell case, we can show F(t eo )/ F^ max j < Tt^j/teo ~ 0.6 
where F(t eo ) is the the optical flux from the reverse shock at the earliest observations by 
Fox (2002) at t eo ~ 0.01 days. This limit is inconsistent with the observations, and hence 

Vm,r(t x ) < Vr. 

If v m ,r(tx) < vr, the reverse shock emission simply decreases as t~ 2 . Since the extrapo- 
lation of the light curve from the data points by Fox (2002) to the earlier time T ~ 10~ 3 day 
with a power law of t~ 2 violates the upper- limits by Torii et al (2002) (see figl), GRB 021004 
should be a thin shell case, in which the peak of the reverse shock emission is delayed, i.e. 
t x >T. 
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Using eqs. (5) and (6) with p = 2.4, t m j ~ 0.06 days, F u>max j ~ 1.3 mJy (R-band), 
E = 5.6 x 10 52 ergs, and z = 2.3, we search for a set of parameters (e B ,rj) with which the 
theoretical estimates give a reasonable fit to all observations. With a normalized Lorentz 
factor k = i]/r] c , we can show that the crossing time t x oc /t~ 8 / 3 and the optical flux from the 
reverse shock at the crossing time F(t x ) = [y B j 1 v m ,r)^' yP ^^' 2 'F v ^ max ^ oc ^ v t~ B p ^ 8 . Since the 
optical flux F(t x ) depends on e B very weakly oc e B 1 ^ 20 (p = 2.4) and e^ 40 (j» = 2.2), e B is not 
well determined from the light curve, or equivalently, we can explain the peculiar behavior 
of the light curve with a wide range of e B . When changing the value of k, the position of the 
peak moves along a line of F oc f 3p / 8 on the flux - time plane, and k ~ 0.55 gives the best 
fit to the observations. When we fit the observations at > 0.1 day with a flatter power law 
of t 3 ( 1_p )/ 4 = t~ 0,9 (p = 2.2), we need to chose a earlier peak time t m j. Assuming the same 
parameters with the p = 2.4 case, except t m j ~ 0.035 days, the best fit is given by k ~ 0.65. 

Since e e and n are determined by eqs. (5) and (6), v m j and F v , max j do not depend on 
e B . The value of e B is constrained only by the the cooling break. The cooling frequency u c j oc 

rV* crosses the optical band at t ~ 0.79 (^) 3 (ff ) 3 fe) 4 ' (^y 2 (ff 

days. Holland et al. (2002a) reported that there is no evidence for color evolution between 
8.5 hours and 5.5 days after the burst (see, however, Matheson et al. 2002; Bersier et al. 
2002), and hence e# < 0.14. A lower limit e B > 8 x 10~ 5 is required from e e < 1. 

Fig 1 shows an example case in which p = 2.4 and e B = 3.0 x 10~ 3 are assumed. This 
choice leads to e e ~ 0.28, n ~ 0.45 proton cm' 3 and rj ~ 120. These are surprisingly 
typical values obtained in other afterglow observations (Panaitescu and Kumar 2002). The 
dashed and dashed dotted lines show the optical light curve of the reverse and forward 
shock emission, respectively. The thick solid lines depicts the total flux. Around the peak 
time of the forward shock ~ 0.1 day, our estimate slightly deviates from the observations. 
However, in our estimate, we assumed a simplified synchrotron spectrum which is described 
by a broken power law. Since a more realistic synchrotron spectrum is rounded at the break 
frequencies (Granot, Piran & Sari 1999), a light curve should be also rounded at the break 
time (dotted line). The short time scale variability, which is prominent around ~ 1 day, 
might be produced by ISM turbulence (Wang & Loeb 2000; Holland et al. 2002b; Lazzati 
et al 2002). The latest data point in Fig 1 is lower than the extrapolation with a scaling of 
£-i.05^ fj^jg mig]-^ be a signature of the jet break. 

Assuming p = 2.2 and e B = 3.0 x 10~ 3 , we get e e ~ 0.32, n ~ 0.45 proton cm -3 and 
i] ~ 140. The total flux is shown by the thin solid line. In this case, we need to assume 
larger ISM turbulence. 
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3. X-ray and Radio Afterglow 

In this section, we assume the values of parameters with which we have shown the 
example cases in the previous section, and estimate the X-ray and radio afterglows. The 
extinction correction a u is unity for X-ray and radio afterglows, and hence F v>max j ~ 1.6 
mJy (radio and X-ray). 

X-ray Afterglow: Since the X-ray band ~ 5keV is well above the typical frequency 
of the reverse shock emission, the contribution from the reverse shock to the X-ray band is 
negligible. The X-ray afterglow should be described only by the forward shock emission. The 
luminosity in X-ray band should decrease as t^ 3p ^ 4 ~ t~ L3 (p = 2.4) or t~ 1A5 (p = 2.2). The 
Chandra X-ray observatory observed the afterglow for a total exposure of 87ksec, beginning 
at Oct 5 8:55 UT (Sako and Harrison 2002). The count rate decrease with a power law slope 
of —1.0 ± 0.2. The mean 2-10 keV X-ray flux is ~ 4.3 x 10~ 13 ergs cm" 2 s _1 . We estimate 
the 5 keV flux at the observational mean time 1.36 days ~ 3.1 x 10~ 13 ergs cm -2 s" 1 for 
p = 2.4 and ~ 6.4 x 10~ 13 ergs cm -2 s _1 for p = 2.2. Our estimates are in a good agreement 
with the observations. 

Radio Afterglow: The forward shock emission in radio band ~ 10 GHz increases as 
t 1 / 2 until the flux reaches to the maximum ~ 1.6 mJy at ~ 80 days for p = 2.4 (dashed 
dotted line in fig2) or at ~ 50 days for p = 2.2. After the typical frequency v m , r crosses the 
radio band, the reverse shock emission decays as ~ t~ 2 (dashed line for p = 2.4). At low 
frequencies and early times, self absorption takes an important role and significantly reduces 
the flux. A simple estimate of the maximal flux (dotted line for p = 2.4) is the emission 
from the black body with the reverse shock temperature (Kobayashi & Sari 2000). The thick 
and thin solid line depicts the total flux for p = 2.4 and for p = 2.2, respectively. Since the 
observations (circles) are done in various frequencies, we scaled the observed value to the 
expected value at 10GHz by using a spectral slope of 1 1 . This burst might also cause a 
bright radio flare ~ 1 mJy around ~ 0.5 day as observed in GRB 990123. When we fit the 
optical observations with p = 2.2, an earlier optical peak time t m j is required. Since the 
peak time at the radio band is proportional to t m j, the modeling with p = 2.2 gives a better 
fit to the radio observations. 



1 Bcrger ct al. (2002) reported an unusual spectral spectral slope F v <~ v between 8.5 GHz and 86GHz 
from the observations with the VLA on October 10.17 UT, and claimed that the spectrum is not due to a 
transition from optically-thick [y 2 ) to optically thin (i/ 1 / 3 ) emission. The superposition of the forward and 
reverse shock emission could give even flatter spectrum. Non-standard emission mechanism may be needed 
to explain this unusual spectrum. 
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4. Conclusions 

The low-frequency (e.g. optical, IR & radio) lightcurve from the forward shock is ex- 
pected to rise initially and to decay after the typical synchrotron frequency crosses the 
observational band. Although such a behavior has been observed in the radio band, pre- 
vious optical afterglow observations were only made at too late times to catch the rising 
phase. However, the swift localization of GRB 021004 by HETE II allowed the follow-up of 
the afterglow at very early time. This burst has so far the earliest detected optical afterglow. 
We have shown that the rising phase of the optical emission might be caught for the first 
time in GRB 021004, together with a reverse shock emission. The superposition of both 
the forward shock and the reverse shock emissions can well fit for the 0.1 day re-brightening 
feature. With the standard values of parameters inferred from other afterglow observations, 
we have constructed example cases in which theoretical estimates fit the broadband obser- 
vations. We therefore suggest that the early re-brightening might be a common feature in 
optical afterglows. 

The reverse shock emission in GRB 990123 peaked at about ~ 1 Jy in R-band, while 
the peak in GRB 021004 was only a few mJy Since GRB 990123 is a very bright burst with 
a fluence about 100 times that of GRB 021004, we expect GRB 990123 could produce the 
much brighter reverse shock emission. Another difference is the typical frequency v m>r (t x ). 
In GRB 990123 it is close to the R-band (Sari & Piran 199b), while in GRB 021004 it is 
estimated as ~ 1.4 x 10 12 Hz (p = 2.4) or ~ 8.8 x 10 11 Hz (p = 2.2). The lower v m , r (t x ) also 
makes the reverse shock emission dimmer in the R-band (Kobayashi 2000). 

We thank Peter Meszaros for valuable comments and the anonymous referee for valuable 
suggestions. We acknowledge support through the Center for Gravitational Wave Physics, 
which is funded by NSF under cooperative agreement PHY 01-14375, and through NASA 
NAG5-9192. 
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Fig. 1. — Optical light curve: (1) modeling with p = 2.4, forward shock emission (dashed 
dotted), reverse shock emission (dashed), total flux (thick solid), accurate spectrum mode 
(dotted). e B = 3.0 x 10~ 3 , e e = 0.28. n = 0.45 protons cm" 3 , r] = 120 and E = 5.6 x 10 52 
ergs. (2) modeling with p = 2.2, total flux (thin solid). e B = 3.0 x 10~ 3 , e e = 0.32. n = 0.45 
protons cm -3 , 77 = 140 and E = 5.6 x 10 52 ergs. Measurements (circles) and upper-limits 
(triangles). Data from: Bersier et al. 2002b; Fox 2002; Halpern et al. 2002a,b; Holland et al. 
2002c,d; Malesani et al. 2002a,b; Masetti et al. 2002; Matsumoto et al. 2002a,b; Mirabal et 
al. 2002a,b; Oksanen&Aho 2002; Oksanen et al. 2002; Sahu et al. 2002; Stanek et al. 2002; 
Stefanon et al. 2002; Torii et al. 2002; Uemura et al. 1566; Weidinger et al. 2002; Winn et 
al. 2002; Zharikov et al. 2002, following the calibration of Henden 2002. 
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Fig. 2. — Radio light curve: (1) modeling with p = 2.4, forward shock emission (dashed 
dotted), reverse shock emission (dashed), self-absorption limit (dotted), total flux (thick 
solid). (2) modeling with p = 2.2, total flux (thin solid). The parameters are the same 
as in Fig.l. Measurements with error bars (circles), measurement without error bar (star) 
and upper-limit (triangle). Data from: Frail & Berger 2002; Pooley 2002a,b,c; Berger et al. 
2002b. 



